Identification of master regulator genes controlling pathogenic CD4+ T cell fate in inflammatory bowel disease through transcriptional network analysis.

Inflammatory bowel diseases (IBD) are a group of chronic inflammatory conditions of the gastrointestinal tract associated with multiple pathogenic factors, including dysregulation of the immune response. Effector CD4+ T cells and regulatory CD4+ T cells (Treg) are central players in maintaining the balance between tolerance and inflammation. Interestingly, genetic modifications in these cells have been implicated in regulating the commitment of specific phenotypes and immune functions. However, the transcriptional program controlling the pathogenic behavior of T helper cells in IBD progression is still unknown. In this study, we aimed to find master transcription regulators controlling the pathogenic behavior of effector CD4+ T cells upon gut inflammation. To achieve this goal, we used an animal model of IBD induced by the transfer of naïve CD4+ T cells into recombination-activating gene 1 (Rag1) deficient mice, which are devoid of lymphocytes. As a control, a group of Rag1-/- mice received the transfer of the whole CD4+ T cells population, which includes both effector T cells and Treg. When gut inflammation progressed, we isolated CD4+ T cells from the colonic lamina propria and spleen tissue, and performed bulk RNA-seq. We identified differentially up- and down-regulated genes by comparing samples from both experimental groups. We found 532 differentially expressed genes (DEGs) in the colon and 30 DEGs in the spleen, mostly related to Th1 response, leukocyte migration, and response to cytokines in lamina propria T-cells. We integrated these data into Gene Regulatory Networks to identify Master Regulators, identifying four up-regulated master gene regulators (Lef1, Dnmt1, Mybl2, and Jup) and only one down-regulated master regulator (Foxo3). The altered expression of master regulators observed in the transcriptomic analysis was confirmed by qRT-PCR analysis and found an up-regulation of Lef1 and Mybl2, but without differences on Dnmt1, Jup, and Foxo3. These two master regulators have been involved in T cells function and cell cycle progression, respectively. We identified two master regulator genes associated with the pathogenic behavior of effector CD4+ T cells in an animal model of IBD. These findings provide two new potential molecular targets for treating IBD.

Pangenome-level analysis of nucleoid-associated proteins in the Acidithiobacillia class: insights into their functional roles in mobile genetic elements biology.

Mobile genetic elements (MGEs) are relevant agents in bacterial adaptation and evolutionary diversification. Stable appropriation of these DNA elements depends on host factors, among which are the nucleoid-associated proteins (NAPs). NAPs are highly abundant proteins that bind and bend DNA, altering its topology and folding, thus affecting all known cellular DNA processes from replication to expression. Even though NAP coding genes are found in most prokaryotic genomes, their functions in host chromosome biology and xenogeneic silencing are only known for a few NAP families. Less is known about the occurrence, abundance, and roles of MGE-encoded NAPs in foreign elements establishment and mobility. In this study, we used a combination of comparative genomics and phylogenetic strategies to gain insights into the diversity, distribution, and functional roles of NAPs within the class Acidithiobacillia with a special focus on their role in MGE biology. Acidithiobacillia class members are aerobic, chemolithoautotrophic, acidophilic sulfur-oxidizers, encompassing substantial genotypic diversity attributable to MGEs. Our search for NAP protein families (PFs) in more than 90 genomes of the different species that conform the class, revealed the presence of 1,197 proteins pertaining to 12 different NAP families, with differential occurrence and conservation across species. Pangenome-level analysis revealed 6 core NAP PFs that were highly conserved across the class, some of which also existed as variant forms of scattered occurrence, in addition to NAPs of taxa-restricted distribution. Core NAPs identified are reckoned as essential based on the conservation of genomic context and phylogenetic signals. In turn, various highly diversified NAPs pertaining to the flexible gene complement of the class, were found to be encoded in known plasmids or, larger integrated MGEs or, present in genomic loci associated with MGE-hallmark genes, pointing to their role in the stabilization/maintenance of these elements in strains and species with larger genomes. Both core and flexible NAPs identified proved valuable as markers, the former accurately recapitulating the phylogeny of the class, and the later, as seed in the bioinformatic identification of novel episomal and integrated mobile elements.

Protein-Mediated Electroporation in a Cardiac Voltage-Sensing Domain Due to an nsPEF Stimulus.

This study takes a step in understanding the physiological implications of the nanosecond pulsed electric field (nsPEF) by integrating molecular dynamics simulations and machine learning techniques. nsPEF, a state-of-the-art technology, uses high-voltage electric field pulses with a nanosecond duration to modulate cellular activity. This investigation reveals a relatively new and underexplored phenomenon: protein-mediated electroporation. Our research focused on the voltage-sensing domain (VSD) of the NaV1.5 sodium cardiac channel in response to nsPEF stimulation. We scrutinized the VSD structures that form pores and thereby contribute to the physical chemistry that governs the defibrillation effect of nsPEF. To do so, we conducted a comprehensive analysis involving the clustering of 142 replicas simulated for 50 ns under nsPEF stimuli. We subsequently pinpointed the representative structures of each cluster and computed the free energy between them. We find that the selected VSD of NaV1.5 forms pores under nsPEF stimulation, but in a way that significant differs from the traditional VSD opening. This study not only extends our understanding of nsPEF and its interaction with protein channels but also adds a new effect to further study.

Nanosecond Pulsed Electric Field (nsPEF): Opening the Biotechnological Pandora's Box.

Nanosecond Pulsed Electric Field (nsPEF) is an electrostimulation technique first developed in 1995; nsPEF requires the delivery of a series of pulses of high electric fields in the order of nanoseconds into biological tissues or cells. They primary effects in cells is the formation of membrane nanopores and the activation of ionic channels, leading to an incremental increase in cytoplasmic Ca2+ concentration, which triggers a signaling cascade producing a variety of effects: from apoptosis up to cell differentiation and proliferation. Further, nsPEF may affect organelles, making nsPEF a unique tool to manipulate and study cells. This technique is exploited in a broad spectrum of applications, such as: sterilization in the food industry, seed germination, anti-parasitic effects, wound healing, increased immune response, activation of neurons and myocites, cell proliferation, cellular phenotype manipulation, modulation of gene expression, and as a novel cancer treatment. This review thoroughly explores both nsPEF's history and applications, with emphasis on the cellular effects from a biophysics perspective, highlighting the role of ionic channels as a mechanistic driver of the increase in cytoplasmic Ca2+ concentration.

Expression of KID syndromic mutation Cx26S17F produces hyperactive hemichannels in supporting cells of the organ of Corti.

Some mutations in gap junction protein Connexin 26 (Cx26) lead to syndromic deafness, where hearing impairment is associated with skin disease, like in Keratitis Ichthyosis Deafness (KID) syndrome. This condition has been linked to hyperactivity of connexin hemichannels but this has never been demonstrated in cochlear tissue. Moreover, some KID mutants, like Cx26S17F, form hyperactive HCs only when co-expressed with other wild-type connexins. In this work, we evaluated the functional consequences of expressing a KID syndromic mutation, Cx26S17F, in the transgenic mouse cochlea and whether co-expression of Cx26S17F and Cx30 leads to the formation of hyperactive HCs. Indeed, we found that cochlear explants from a constitutive knock-in Cx26S17F mouse or conditional in vitro cochlear expression of Cx26S17F produces hyperactive HCs in supporting cells of the organ of Corti. These conditions also produce loss of hair cells stereocilia. In supporting cells, we found high co-localization between Cx26S17F and Cx30. The functional properties of HCs formed in cells co-expressing Cx26S17F and Cx30 were also studied in oocytes and HeLa cells. Under the recording conditions used in this study Cx26S17F did not form functional HCs and GJCs, but cells co-expressing Cx26S17F and Cx30 present hyperactive HCs insensitive to HCs blockers, Ca2+ and La3+, resulting in more Ca2+ influx and cellular damage. Molecular dynamic analysis of putative heteromeric HC formed by Cx26S17F and Cx30 presents alterations in extracellular Ca2+ binding sites. These results support that in KID syndrome, hyperactive HCs are formed by the interaction between Cx26S17F and Cx30 in supporting cells probably causing damage to hair cells associated to deafness.

On the effectiveness of communication strategies as non-pharmaceutical interventions to tackle epidemics.

When pharmaceutical interventions are unavailable to deal with an epidemic outbreak, adequate management of communication strategies can be key to reduce the contagion risks. On the one hand, accessibility to trustworthy and timely information, whilst on the other, the adoption of preventive behaviors may be both crucial. However, despite the abundance of communication strategies, their effectiveness has been scarcely evaluated or merely circumscribed to the scrutiny of public affairs. To study the influence of communication strategies on the spreading dynamics of an infectious disease, we implemented a susceptible-exposed-infected-removed-dead (SEIRD) epidemiological model, using an agent-based approach. Agents in our systems can obtain information modulating their behavior from two sources: (i) through the local interaction with other neighboring agents and, (ii) from a central entity delivering information with a certain periodicity. In doing so, we highlight how global information delivered from a central entity can reduce the impact of an infectious disease and how informing even a small fraction of the population has a remarkable impact, when compared to not informing the population at all. Moreover, having a scheme of delivering daily messages makes a stark difference on the reduction of cases, compared to the other evaluated strategies, denoting that daily delivery of information produces the largest decrease in the number of cases. Furthermore, when the information spreading relies only on local interactions between agents, and no central entity takes actions along the dynamics, then the epidemic spreading is virtually independent of the initial amount of informed agents. On top of that, we found that local communication plays an important role in an intermediate regime where information coming from a central entity is scarce. As a whole, our results highlight the importance of proper communication strategies, both accurate and daily, to tackle epidemic outbreaks.

Exploring the Conformational Changes Induced by Nanosecond Pulsed Electric Fields on the Voltage Sensing Domain of a Ca2+ Channel.

Nanosecond Pulsed Electric Field (nsPEF or Nano Pulsed Stimulation, NPS) is a technology that delivers a series of pulses of high-voltage electric fields during a short period of time, in the order of nanoseconds. The main consequence of nsPEF upon cells is the formation of nanopores, which is followed by the gating of ionic channels. Literature is conclusive in that the physiological mechanisms governing ion channel gating occur in the order of milliseconds. Hence, understanding how these channels can be activated by a nsPEF would be an important step in order to conciliate fundamental biophysical knowledge with improved nsPEF applications. To get insights on both the kinetics and thermodynamics of ion channel gating induced by nsPEF, in this work, we simulated the Voltage Sensing Domain (VSD) of a voltage-gated Ca2+ channel, inserted in phospholipidic membranes with different concentrations of cholesterol. We studied the conformational changes of the VSD under a nsPEF mimicked by the application of a continuous electric field lasting 50 ns with different intensities as an approach to reveal novel mechanisms leading to ion channel gating in such short timescales. Our results show that using a membrane with high cholesterol content, under an nsPEF of 50 ns and E→ = 0.2 V/nm, the VSD undergoes major conformational changes. As a whole, our work supports the notion that membrane composition may act as an allosteric regulator, specifically cholesterol content, which is fundamental for the response of the VSD to an external electric field. Moreover, changes on the VSD structure suggest that the gating of voltage-gated Ca2+ channels by a nsPEF may be due to major conformational changes elicited in response to the external electric field. Finally, the VSD/cholesterol-bilayer under an nsPEF of 50 ns and E→ = 0.2 V/nm elicits a pore formation across the VSD suggesting a new non-reported effect of nsPEF into cells, which can be called a "protein mediated electroporation".

Simulations on Simple Models of Connexin Hemichannels Indicate That Ca2+ Blocking Is Not a Pure Electrostatic Effect.

Connexin hemichannels allow the unspecific but regulated interchange of molecules from ions to second messenger and ATP, between the eukariotic cell and its extracellular space. The transport of ions and water through hemichannels is important for physiological functions and also in the progression of several pathological conditions. Extracellular Ca2+ concentration is one of the regulators that drives the channel to a closed state. However the relation between their functional and structural states is far for being totally understood. In this work, we modelled connexin hemichannels using simple systems based on a fixed array of carbon atoms and assess the Ca2+ regulation using molecular dynamics simulations. The two proposed mechanism described so far for calcium action were studied combined, e.g., an electrostatic effect and a pore stretching. Our results show that the addition of positive charge density inside the channel cannot stop the flow of potassium, chloride nor water. Only a pore stretching at the center of the pore can explain the channel blocking.

Brief Comparison between Experimental and Computationally Generated Ensembles of RNA Dinucleotides.

The computational modeling of RNA and its interactions is of crucial importance for the understanding of the wide variety of biological functions it performs. Among these approaches, several coarse-grained models employ structural databases to derive their energy functions or to define scoring functions for structure prediction purposes. In many cases, the parametrization is done by using as a reference a set of experimentally determined structures or data obtained from Molecular Dynamics simulations. Since the two choices are clearly different, we present here a brief comparison of the essential spaces of a set of conformations of two topologically connected nucleotides generated by these means. We find that when the nucleotides are embedded into a duplex, the essential spaces of both ensembles are quite restricted and do not differ substantially. Nevertheless, when the conformational space of a free dinucleoside monophosphate simulation is compared against a similar set obtained from the structural database, the differences of the essential spaces are considerable. In addition, we show a brief comparison of a specific distance between the nucleotides which correlates with the sugar pucker of the first nucleotide and analyze its distribution under similar conditions.

l-DOPA modulates the kinetics but not the thermodynamic equilibrium of TTA+ amphiphiles forming lyotropic nematic liquid crystals.

Lyotropic liquid crystals (LLCs) are mixtures of amphiphile molecules usually studied as mimetic of biological membrane. The equilibrium dynamics of tetradecyltrimethyl ammonium cation (TTA+) molecules forming nematic LLCs (LNLCs) is guided by a dive-in mechanism where TTA+ molecules spontaneously leave and re-enter the bicelle. Of note, this dynamic behavior could be exploited to produce drug nano-delivery systems based on LNLCs. Therefore, the understanding of the effect of pharmaceutically interesting molecules in the dynamics of the dive-in mechanism should be crucial for drug delivery applications. In this work, we studied the effects of l-DOPA in the equilibrium dynamics of TTA+ bicelles forming LNLCs, employing a transdisciplinary approach based on 2H-NMR together with molecular modeling and molecular dynamics simulations. Our data suggest that l-DOPA perturbs the kinetic of the dive-in mechanism but not the thermodynamics of this process. As whole, our results provide fundamental insights on the mechanisms by which l-DOPA govern the equilibrium of LNLCs bicelles.

Orientational and Folding Thermodynamics via Electric Dipole Moment Restraining.

The projection of molecular processes onto a small set of relevant descriptors, the so-called reaction coordinates or collective variables (CVs), is a technique nowadays routinely employed by the biomolecular simulation community. In this work, we implemented two CVs to manipulate the orientation (i.e., angle) (µâƒ—a) and magnitude (|µâƒ—|) of the electric dipole moment. In doing so, we studied the thermodynamics of water orientation under the application of external voltages and the folding of two polypeptides at zero-field conditions. The projection of the free-energy [potential of mean force (PMF)] along water orientation defined an upper limit of around 0.3 V for irrelevant thermodynamic effects. On the other hand, sufficiently strong µâƒ—a restraints applied on 12-alanine (Ala12) triggered structural effects because of the alignment of local dipoles; for lower restraints, a full-body rotation is achieved. The manipulation of |µâƒ—| produced strong perturbations on the secondary structure of Ala12, promoting an enhanced sampling to its configurational space. Rigorous free-energy calculations in the form of 2-D PMFs for deca-alanine showed the utility of |µâƒ—| as a reaction coordinate to study folding in small α helices. As a whole, we propose that the manipulation of both components of the dipole moment, µâƒ—a and |µâƒ—|, provides thermodynamics insights into the structural conformation and stability of biomolecules. These new CVs are implemented in the Colvars module, available for NAMD and LAMMPS.

RIP-MD: a tool to study residue interaction networks in protein molecular dynamics.

Protein structure is not static; residues undergo conformational rearrangements and, in doing so, create, stabilize or break non-covalent interactions. Molecular dynamics (MD) is a technique used to simulate these movements with atomic resolution. However, given the data-intensive nature of the technique, gathering relevant information from MD simulations is a complex and time consuming process requiring several computational tools to perform these analyses. Among different approaches, the study of residue interaction networks (RINs) has proven to facilitate the study of protein structures. In a RIN, nodes represent amino-acid residues and the connections between them depict non-covalent interactions. Here, we describe residue interaction networks in protein molecular dynamics (RIP-MD), a visual molecular dynamics (VMD) plugin to facilitate the study of RINs using trajectories obtained from MD simulations of proteins. Our software generates RINs from MD trajectory files. The non-covalent interactions defined by RIP-MD include H-bonds, salt bridges, VdWs, cation-π, π-π, Arginine-Arginine, and Coulomb interactions. In addition, RIP-MD also computes interactions based on distances between Cαs and disulfide bridges. The results of the analysis are shown in an user friendly interface. Moreover, the user can take advantage of the VMD visualization capacities, whereby through some effortless steps, it is possible to select and visualize interactions described for a single, several or all residues in a MD trajectory. Network and descriptive table files are also generated, allowing their further study in other specialized platforms. Our method was written in python in a parallelized fashion. This characteristic allows the analysis of large systems impossible to handle otherwise. RIP-MD is available at http://www.dlab.cl/ripmd.

Rule-Based Models and Applications in Biology.

Complex systems are governed by dynamic processes whose underlying causal rules are difficult to unravel. However, chemical reactions, molecular interactions, and many other complex systems can be usually represented as concentrations or quantities that vary over time, which provides a framework to study these dynamic relationships. An increasing number of tools use these quantifications to simulate dynamically complex systems to better understand their underlying processes. The application of such methods covers several research areas from biology and chemistry to ecology and even social sciences.In the following chapter, we introduce the concept of rule-based simulations based on the Stochastic Simulation Algorithm (SSA) as well as other mathematical methods such as Ordinary Differential Equations (ODE) models to describe agent-based systems. Besides, we describe the mathematical framework behind Kappa (κ), a rule-based language for the modeling of complex systems, and some extensions for spaßtial models implemented in PISKaS (Parallel Implementation of a Spatial Kappa Simulator). To facilitate the understanding of these methods, we include examples of how these models can be used to describe population dynamics in a simple predator-prey ecosystem or to simulate circadian rhythm changes.

The syndromic deafness mutation G12R impairs fast and slow gating in Cx26 hemichannels.

Mutations in connexin 26 (Cx26) hemichannels can lead to syndromic deafness that affects the cochlea and skin. These mutations lead to gain-of-function hemichannel phenotypes by unknown molecular mechanisms. In this study, we investigate the biophysical properties of the syndromic mutant Cx26G12R (G12R). Unlike wild-type Cx26, G12R macroscopic hemichannel currents do not saturate upon depolarization, and deactivation is faster during hyperpolarization, suggesting that these channels have impaired fast and slow gating. Single G12R hemichannels show a large increase in open probability, and transitions to the subconductance state are rare and short-lived, demonstrating an inoperative fast gating mechanism. Molecular dynamics simulations indicate that G12R causes a displacement of the N terminus toward the cytoplasm, favoring an interaction between R12 in the N terminus and R99 in the intracellular loop. Disruption of this interaction recovers the fast and slow voltage-dependent gating mechanisms. These results suggest that the mechanisms of fast and slow gating in connexin hemichannels are coupled and provide a molecular mechanism for the gain-of-function phenotype displayed by the syndromic G12R mutation.

Stochastic simulation of multiscale complex systems with PISKaS: A rule-based approach.

Computational simulation is a widely employed methodology to study the dynamic behavior of complex systems. Although common approaches are based either on ordinary differential equations or stochastic differential equations, these techniques make several assumptions which, when it comes to biological processes, could often lead to unrealistic models. Among others, model approaches based on differential equations entangle kinetics and causality, failing when complexity increases, separating knowledge from models, and assuming that the average behavior of the population encompasses any individual deviation. To overcome these limitations, simulations based on the Stochastic Simulation Algorithm (SSA) appear as a suitable approach to model complex biological systems. In this work, we review three different models executed in PISKaS: a rule-based framework to produce multiscale stochastic simulations of complex systems. These models span multiple time and spatial scales ranging from gene regulation up to Game Theory. In the first example, we describe a model of the core regulatory network of gene expression in Escherichia coli highlighting the continuous model improvement capacities of PISKaS. The second example describes a hypothetical outbreak of the Ebola virus occurring in a compartmentalized environment resembling cities and highways. Finally, in the last example, we illustrate a stochastic model for the prisoner's dilemma; a common approach from social sciences describing complex interactions involving trust within human populations. As whole, these models demonstrate the capabilities of PISKaS providing fertile scenarios where to explore the dynamics of complex systems.

LoTo: a graphlet based method for the comparison of local topology between gene regulatory networks.

One of the main challenges of the post-genomic era is the understanding of how gene expression is controlled. Changes in gene expression lay behind diverse biological phenomena such as development, disease and the adaptation to different environmental conditions. Despite the availability of well-established methods to identify these changes, tools to discern how gene regulation is orchestrated are still required. The regulation of gene expression is usually depicted as a Gene Regulatory Network (GRN) where changes in the network structure (i.e., network topology) represent adjustments of gene regulation. Like other networks, GRNs are composed of basic building blocks; small induced subgraphs called graphlets. Here we present LoTo, a novel method that using Graphlet Based Metrics (GBMs) identifies topological variations between different states of a GRN. Under our approach, different states of a GRN are analyzed to determine the types of graphlet formed by all triplets of nodes in the network. Subsequently, graphlets occurring in a state of the network are compared to those formed by the same three nodes in another version of the network. Once the comparisons are performed, LoTo applies metrics from binary classification problems calculated on the existence and absence of graphlets to assess the topological similarity between both network states. Experiments performed on randomized networks demonstrate that GBMs are more sensitive to topological variation than the same metrics calculated on single edges. Additional comparisons with other common metrics demonstrate that our GBMs are capable to identify nodes whose local topology changes between different states of the network. Notably, due to the explicit use of graphlets, LoTo captures topological variations that are disregarded by other approaches. LoTo is freely available as an online web server at http://dlab.cl/loto.

On Biophysical Properties and Sensitivity to Gap Junction Blockers of Connexin 39 Hemichannels Expressed in HeLa Cells.

Although connexins (Cxs) are broadly expressed by cells of mammalian organisms, Cx39 has a very restricted pattern of expression and the biophysical properties of Cx39-based channels [hemichannels (HCs) and gap junction channels (GJCs)] remain largely unknown. Here, we used HeLa cells transfected with Cx39 (HeLa-Cx39 cells) in which intercellular electrical coupling was not detected, indicating the absence of GJCs. However, functional HCs were found on the surface of cells exposed to conditions known to increase the open probability of other Cx HCs (e.g., extracellular divalent cationic-free solution (DCFS), extracellular alkaline pH, mechanical stimulus and depolarization to positive membrane potentials). Cx39 HCs were blocked by some traditional Cx HC blockers, but not by others or a pannexin1 channel blocker. HeLa-Cx39 cells showed similar resting membrane potentials (RMPs) to those of parental cells, and exposure to DCFS reduced RMPs in Cx39 transfectants, but not in parental cells. Under these conditions, unitary events of ~75 pS were frequent in HeLa-Cx39 cells and absent in parental cells. Real-time cellular uptake experiments of dyes with different physicochemical features, as well as the application of a machine-learning approach revealed that Cx39 HCs are preferentially permeable to molecules characterized by six categories of descriptors, namely: (1) electronegativity, (2) ionization potential, (3) polarizability, (4) size and geometry, (5) topological flexibility and (6) valence. However, Cx39 HCs opened by mechanical stimulation or alkaline pH were impermeable to Ca2+. Molecular modeling of Cx39-based channels suggest that a constriction present at the intracellular portion of the para helix region co-localizes with an electronegative patch, imposing an energetic and steric barrier, which in the case of GJCs may hinder channel function. Results reported here demonstrate that Cx39 form HCs and add to our understanding of the functional roles of Cx39 HCs under physiological and pathological conditions in cells that express them.

Graphlet Based Metrics for the Comparison of Gene Regulatory Networks.

Understanding the control of gene expression remains one of the main challenges in the post-genomic era. Accordingly, a plethora of methods exists to identify variations in gene expression levels. These variations underlay almost all relevant biological phenomena, including disease and adaptation to environmental conditions. However, computational tools to identify how regulation changes are scarce. Regulation of gene expression is usually depicted in the form of a gene regulatory network (GRN). Structural changes in a GRN over time and conditions represent variations in the regulation of gene expression. Like other biological networks, GRNs are composed of basic building blocks called graphlets. As a consequence, two new metrics based on graphlets are proposed in this work: REConstruction Rate (REC) and REC Graphlet Degree (RGD). REC determines the rate of graphlet similarity between different states of a network and RGD identifies the subset of nodes with the highest topological variation. In other words, RGD discerns how th GRN was rewired. REC and RGD were used to compare the local structure of nodes in condition-specific GRNs obtained from gene expression data of Escherichia coli, forming biofilms and cultured in suspension. According to our results, most of the network local structure remains unaltered in the two compared conditions. Nevertheless, changes reported by RGD necessarily imply that a different cohort of regulators (i.e. transcription factors (TFs)) appear on the scene, shedding light on how the regulation of gene expression occurs when E. coli transits from suspension to biofilm. Consequently, we propose that both metrics REC and RGD should be adopted as a quantitative approach to conduct differential analyses of GRNs. A tool that implements both metrics is available as an on-line web server (http://dlab.cl/loto).

Exploring the Membrane Potential of Simple Dual-Membrane Systems as Models for Gap-Junction Channels.

The conductance of ion channels can be modulated by a transmembrane potential difference, due to alterations on ion-mobility and also by changes in the pore structure. Despite the vast knowledge regarding the influence of voltage on transport properties of ion channels, little attention has been paid to describe, with atomic detail, the modulation of ionic transport in gap-junction channels (GJCs). Hence, molecular dynamics simulations were performed to explore the conductance of simple dual-membrane systems that account for the very basic features of GJCs. In doing so, we studied the influence of different charge distributions in the channel surface on these idealized systems under external electric fields, paying attention to the behavior of the electrostatic potential, ion density, ion currents, and equilibrium properties. Our results demonstrate that the incorporation of a charge distribution akin GJCs decreased anionic currents, favoring the transport of cationic species. Moreover, a thermodynamic characterization of ionic transport in these systems demonstrate the existence of a kinetic barrier that hinders anionic currents, reinforcing the role played by the internal arrangement of charges in GJCs. Overall, our results provide insights at the atomic scale on the effects of charge distributions over ionic transport, constituting a step forward into a better understanding of GJCs.

Changes in Olfactory Receptor Expression Are Correlated With Odor Exposure During Early Development in the zebrafish (Danio rerio).

We have previously shown that exposure to phenyl ethyl alcohol (PEA) causes an increase in the expression of the transcription factor otx2 in the olfactory epithelium (OE) of juvenile zebrafish, and this change is correlated with the formation of an odor memory of PEA. Here, we show that the changes in otx2 expression are specific to ßPEA: exposure to αPEA did not affect otx2 expression. We identified 34 olfactory receptors (ORs) representing 16 families on 4 different chromosomes as candidates for direct regulation of OR expression via Otx2. Subsequent in silico analysis uncovered Hnf3b binding sites closely associated with Otx2 binding sites in the regions flanking the ORs. Analysis by quantitative polymerase chain reaction and RNA-seq of OR expression in developing zebrafish exposed to different isoforms of PEA showed that a subset of ORs containing both Otx2/Hnf3b binding sites were downregulated only in ßPEA-exposed juveniles and this change persisted through adult life. Localization of OR expression by in situ hybridization indicates the downregulation occurs at the level of RNA and not the number of cells expressing a given receptor. Finally, analysis of immediate early gene expression in the OE did not reveal changes in c-fos expression in response to either αPEA or ßPEA.